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Abstract 

We study the role of magnetic monopoles at high enough temperature T > 
2Tc, when they can be considered heavy, rare objects embedded into matter 
consisting mostly of the usual "electric" quasiparticles, quarks and gluons. 
We review available lattice results on monopoles at finite temperatures. Then 
we proceed to classical and quantum charge-monopole scattering, solving the 
problem of gluon-monopole scattering for the first time. We find that, while 
this process hardly influences thermodynamic quantities, it does produce a 
large transport cross section, significantly exceeding that for pQCD gluon- 
gluon scattering up to quite high T. Thus, in spite of their relatively small 
density at high T, monopoles are extremely important for QGP transport 
properties, keeping viscosity small enough for hydrodynamics to work at LHC. 

1 Introduction 
1.1 Overview 

As it is known from 1970's, QCD at high temperature T is weakly coupled [1] 
and provides perturbative screening of the charge [2], thus being called Quark- 
Gluon Plasma (QGP). Creating and studying this phase of matter in the laboratory 
has been the goal of experiments at CERN SPS and recently at the Relativistic 
Heavy Ion Collider (RHIC) facility in Brookhaven National Laboratory, soon to 
be continued by the ALICE collaboration at the Large Hadron Collider (LHC). 
RHIC experiments have revealed robust collective phenomena in the form of radial 
and elliptic flows, which turned out to be quite accurately described by near-ideal 
hydrodynamics. QGP thus seems to be the most perfect liquid known, with the 
smallest viscosity-to-entropy ratio rj/s. 



The theory of QGP has shifted from the perturbative-based one, appropriate 
for a weakly coupled (gas) regime, to the non-perturbative methods needed to ad- 
dress the strongly coupled QGP (sQGP for short) regime. This "paradigm shift", 
documented in Refs. P H], is still profoundly affecting the developments. The 
methods addressing strongly coupled gauge theories include in particular the so- 
called AdS/CFT correspondence, relating strongly coupled gauge theory to weakly 
coupled string theory in a particular setting. We will not discuss it in this paper, 
for a recent review see e.g. [Sj. On pure phenomenological grounds, it has been ar- 
gued that, since many substances exhibit a minimum of the viscosity at some phase 
transitions, perhaps QGP is the "best hquid" at the QCD phase transition as well, 
namely at T = Tc [6]. 

Another duality which has been used to explain unusual properties of the sQGP 
is the electric — magnetic duality. Liao and one of us have proposed the so-called 
"magnetic scenario" |7j, according to which the near-Tc region is dominated by 
magnetic monopoles. This is not surprising, if the deconfinement phase transition 
is basically interpreted as their Bose condensation. Based on molecular dynamics of 
classical plasmas with both electric and magnetic quasiparticles, it has been further 
argued in that work that the minimal viscosity /entropy ratio (the "best liquid") does 
not correspond to the phase transition point T = Tc, but rather to the "electric- 
magnetic equilibrium", at T 1.4Tc, where both components of QGP contribute 
about equally to transport coefficients. We will review more recent works on the 
subject below. 

One of the central questions is how sQGP with "perfect fluidity" will change 
into a weakly coupled wQGP with increasing T. In view of the next round of heavy 
ion experiments at LHC, a quite urgent question is what transport properties are 
expected to be observed there, at temperatures reaching about twice those reached 
at RHIC. In order to answer this question, one of course has to understand where 
the "perfect fluidity" property of QGP comes from. As an important example of 
perturbative point of view, we mention the work by Xu, C. Greiner and Stocker [S] 
who argued that the QGP is only moderately coupled, with as = 0.3. .0.6, explaining 
the small viscosity by inclusion of the next-order radiative processes, gg ^ ggg. 
We will discuss this issue partially in the next section, dealing with parametric 
dependences of densities and scattering rates, and also at the end of the paper in 
the discussion section. Here we only notice that, if this should be the explanation, 
one would expect a very slow transition to weakly coupled QGP, induced by the 
logarithmic running of the coupling. 

In this paper we address the issue of QGP transport properties in the "mag- 
netic scenario" framework, moving away from the phase transition region to higher 
temperatures, where QGP is still dominated by the usual electric quasiparticles - 
quarks and gluons - and the coupling is moderately small. Our goal is to study the 
interaction between electric and magnetic sectors. The main difference, with respect 
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to the Liao-Shuryak paper [7], is that their classical (manybody) treatment of this 
interaction (via Lorentz force) is replaced by our (two-body) quantum-mechanical 
calculation. Our main result is the explicit solution of the rather difficult problem 
of quantum gluon-monopole scattering, from which we calculate the corresponding 
transport cross sections. 

Two types of "enhancement mechanisms" are found: one is related to the cou- 
plings and the other one (of much more nontrivial origin) is due to enhanced large- 
angle scattering. The former is rather easy to explain: while "electric" gg, qg, gg 
scattering has cross sections proportional to small dliectric ^ 1; scattering of 
electric-magnetic type contains the product a eiectricO: magnetic which, according to the 
celebrated Dirac condition, cannot run and is equal to an integer (just one for the 
't Hooft-Polyakov monopole we will consider). Therefore, in this problem there is 
no small parameter at any temperature. For many other details the reader should 
read the paper further; nevertheless, we emphasize here our main result, namely 
that the gluon-monopole scattering rates are very large. Parametrically they are 
Wgm/T ~ log(T)~'^, just one power of the log down compared to the gluon-gluon 
pQCD scattering rates Wgg/T ~ log(T)~^. Thus, if the gluon-monopole scattering 
is the dominant one, one should expect stronger decrease of scattering [rj/s growth) 
at LHC, compared to the perturbative scenario 

Below we include a rather extensive introduction to the electric-magnetic dual- 
ity in gauge theories in general, and to the magnetic sector of QCD at high T in 
particular. This choice, which we hope will be useful to many readers, can be justi- 
fied by the somewhat intermittent nature of this field, with many important results 
obtained and then semi-forgotten outside a relatively small group of experts. In par- 
ticular, the scattering of an electric particle on a magnetic monopole is a century-old 
problem. It is hard to explain the quantum scattering without introducing classical 
notations and results first. Thus, although this part of the story is rather old, we 
include some of this material for completeness as an extended introduction. We will 
then proceed to the extension of this theory to gauge particles. 

1.2 Scattering on monopole in classical approximation 

Let us start with a very old (19th century) problem in classical electrodynamics: 
an electrically charged particle with charge e is moving in the magnetic field of a 
static monopole with magnetic coupling g. Classically, one does not need the vector 
potential, thus many subtleties are absent. The interaction is simply given by the 
Lorentz force 

■J, r X f 

mr = —eg — — (1) 

where x indicates the vector product and f = r/r is the unit vector along the line 
connecting both charges. It is worth noticing that only the product of the couplings 
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(eg) appears. 

Furthermore, the cross product of the magnetic field of the monopole and electric 
field of the charge leads to a nonzero Poynting vector, which rotates around r: thus, 
the field itself has a nonzero angular momentum. The total conserved angular 
momentum for this problem has two parts 

J = mf X f + egf. (2) 

The traditional potential scattering only has the first part: therefore, in that case, 
the motion entirely takes place in the so-called "reaction plane" normal to J. In 
the charge-monopole problem, the second term restricts the motion to the so-called 
Poincare cone [H] : its half-opening angle 7r/2 — ^ being 

sin(0 = ^. (3) 

Only at large J (large impact parameter scattering) the angle ^ is small and thus 
the cone opens up, approaching the scattering plane. 

Following Ref. [lO], one can project the motion on the cone to a planar motion, 
by introducing 

R=^—[f-Jif.J)] (4) 
cos(0 

where the first scale factor is introduced to keep the same length for both vectors 
= f^2^ Now, two integrals of motion are 

J = mRx R (5) 
^ mR^ (eqV , ^ 

and the problem seems to be reduced to the motion of a particle of mass m in an 
inverse-square potential. The scattering angle Aip for this planar problem can be 
readily found: it is the variation of ^ as i? goes from oo to its minimum b and back 
to oo 

Note that at large b (small ^) we have Aip ~ as expected for the inverse-square 
potential. Yet this is not the scattering angle of the original problem, because one 
has to project the motion back to the Poincare cone. The true scattering angle - 
namely the angle between the initial and final velocities - is cos 6' = — (£>« ■ Vf). By 
relating velocities on the plane and on the cone one can find it to be 

, .^2( ■ ^ 

cos - = (cos4 j I sm 



2) ^ ^' y 2cos^ 



4 



Thus for distant scattering - small ^ - one gets 6^2^ = 2eg/{mvb), which is much 
larger than Aip ~ 1/6^. The important lesson that we learn from these formulae is 
that the small scattering angle is given by the opening angle of the cone, rather than 
the scattering angle in the planar, inverse- square effective potential. Calculating the 
cross section by da = 27[bdb one finds that, at small angles, it is 

^=^^V- (9) 

dn \mv ) e^' ^ ' 

similar to the Rutherford scattering of two charges. The difference (apart from 
different charges) is also the additional second power of velocity, originating from 
the Lorentz force. 



1.3 Quantum charge-monopole scattering problem 

Jumping from 1900's to 1930's and from classical to quantum mechanics, one first 
has to deal with the necessity of introducing the vector potential A^{x). If this is 
done, the divergence of the magnetic field is zercQ Dirac [TT] brilliantly solved this 
difficulty by introducing the celebrated Dirac string: he then proved it to be an 
unphysical gauge artefact, provided the charge quantization condition is fulfilled 

— = integer. (10) 
47r 

The same condition appears, if we recognize the necessity of quantizing the field 
angular momentum to (semi-integer) multiples of h. 

Now we jump to mid-1970's, when the quantum scattering problem was solved 
in Refs. [ini [E] for the scattering of a scalar particle on a monopole, and in Ref. 
[13] for a monopole-spin 1/2 particle scattering. 

The wave function in spherical coordinates is as usual a sum of products of 
certain r- dependent radial functions, times the angular functions. The former 
basically follow from the inverse- square law potential and thus are easily solved in 
Bessel functions: they correspond to the auxiliary planar projection of the classical 
problem of the previous subsection. The nontrivial part happens to be in the unusual 
angular functions. Before introducing those, let us hint why the usual set of angular 
harmonics Yim{0, should not to be used. The reason is that their classical limit - 
for large values of the indices m ^ / ^ 1 - corresponds to fiat planar motion near 
the z = plane. Nevertheless, we have already learned from the classical limit to 
expect the motion to be concetrated around the Poincare cone instead! 



^ One can introduce the so-called dual vector potential related to B, E in the opposite way: 
then it is the electric field which would be forced to have zero divergence and Dirac strings. 
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The functions that we need (for example in the scalar sector) must satisfy the 
following set of conditions 



f ■ I 




where T is the total angular momentum T = L + /, L is the orbital angular momen- 
tum and / the isospin. The unusual condition in the above set is the last one, since 
the vector J must be projected to the (space-dependent) radial unit vector. The 
functions satisfying this requirement 112] will be introduced in Sec.|3]in the case 
of a scalar particle, and in Sec. |4]in the case of a vector one. Here we can anticipate 
that, in order to check if the hint we provided above is satisfied, it is enough to 
explore the large Z,n limit of the D-functions involved. The result 



exp[-/(0-r)V2] 



(12) 



where cos(6'*) = n/l, shows that they indeed correspond to the Poincare cone. 



1.4 Brief review on QCD monopoles 

The discussion in the previous two subsections was restricted to electrodynamics, in 
which there is neither an explicit explanation of the monopole structure nor a real 
necessity to have themj^ However, non-Abelian gauge theories do have monopole-like 
solitons, as demonstrated by 't Hooft and Polyakov [T^ [T5] in the Georgi-Glashow 
model framework. Three decades later, Seiberg and Witten [T6] have discovered 
beautiful ways in which monopoles and dyons can get alive and even replace the 
usual electric particles (gluons/gluinos) in M=2 Supersymmetric- Yang- Mills theory. 

Since those developments are very well documented in books and reviews, we 
directly jump to finite-T QCD. Explicit calculations show that the static magnetic 
field is not screened perturbatively in QGP [2j: yet Polyakov conjectured |17j that 
magnetic screening should appear non-perturbatively, at the so-called "magnetic 
scale"0EM ~ e^T. Linde |I8] related this magnetic scale and screening to monopole 
solutions, which appear at finite T when the Higgs fundamental scalar of the Georgi- 
Glashow model is replaced by Aq, the component of the gauge field. Shortly 
afterwards, the monopole Bose-condensation was the central idea of the so called 

^The famous Dirac's quote is that he would be surprised that Nature would not make use of 
them. 

In order to keep with the notations used in monopole problems, in this work we use e to 
indicate the usual gauge coupling, reserving g for magnetic coupling. We hope e will not be 
confused with QED charge, which is not used in this work. 
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"color superconductor" picture of confinement, proposed by Mandelstam and 't 
Hooft [in]. Tlie abelian Higgs model (an incarnation of Landau-Ginzburg action) 
was worked out in detail and shown to be in excellent agreement with lattice results 
on the structure of confining flux tubes at T = by a large number of people. The 
profiles of the flux tubes nicely follow the celebrated Abrikosov solution, and lattice 
monopoles do rotate around them, providing a "coil" needed to contain the electric 
flux. 

Lattice monopoles are defined by the procedure [20j which basically locates the 
ends of singular Dirac strings by calculating the total magnetic flux through the 
boundary of elementary 3-d boxes. Since the strings and thus the procedure depend 
on a certain gauge, for decades sceptics kept the viewpoint that those objects are just 
unphysical UV gauge noise. Yet, many specific questions - e.g. monopole density 
and correlations |i21j - produced very reasonable and consistent answers, which are 
apparently independent of the particular lattice parameters. It is hard to imagine 
that lattice gauge artefact have consistent correlations at relatively large distances, 
accurately reproducing those in a Coulomb plasma. We thus assume instead, that 
these monopoles are indeed meaningful physical degrees of freedom, present in QGP 
as quasiparticles and being the source of a Coulomb-like magnetic field. 

The picture of the (uncondensed) monopole gas at T > Tc was discussed in par- 
ticular in the important paper by Smit et al. [22j, in which the relation between 
the monopole density and such observables as the spatial string tension were consid- 
ered. One question (especially important for our work) is whether monopoles form 
a gas or a liquid. To that purpose, the authors of Ref. \12\ have pointed out that 
the combination rtm/M'^ ~ .03 -C 1 is smalQ while in a weakly coupled (Debye) 
plasma it should be large. (Indeed, in high-T electric plasma n ~ T^, Ee ~ eT, 
and thus n/E% ~ 1/e^ ^ 1). Yet such ratio can be small in a liquid; in fact the 
ensemble of instantons in the QCD vacuum - the "instanton liquid" extensively 
studied by one of us [231 [21] - has the same property. The T = instanton density 
^inst ~ 1 fm""^ while the screening mass of the topological charge is m^/ ~ .94 GeV, 
thus riinst/M^, ~ .002 ^ 1 as well. Taking the corresponding roots for comparison, 

1 /3 

one even finds the small parameters in both cases to be similar: rim /Mm ~ .3, 
n]itjM,> ^ .2. 

Recent interest in QGP monopoles (at T > Tc) started with the "magnetic sce- 
nario" [7] which basically suggested the magnetic sector to be a strongly coupled 
(magnetic) Coulomb plasma of monopoles, in its liquid form. Another line of work 
based on lattice monopoles has led Chernodub and Zakharov [25] at the same time 
to a very similar conclusion. An important feature of this scenario [7] is the oppo- 
site running of the electric coupling e and the magnetic one g, induced by the Dirac 
condition eg = const. As recently shown in [26], this feature has been dramatically 

^ For more recent discussion of this feature see Korthals Altes [?f] . 
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confirmed by the behavior of the lattice correlation functions [21] , which indeed dis- 
play monopole-monopole and antimonopole-monopole correlations increasing with 
T. As shown in [26j, those correlations are well described by a picture of classical 
Coulomb gas. The main input in those calculations is the magnetic Coulomb cou- 
pling g growing with T. Furthermore, it was shown to be simply the inverse of the 
gauge coupling e, as Dirac predicted. We consider the correlations observed in [21] 
to be a decisive confirmation of the existence of the long-distance magnetic Coulomb 
field of the monopoles. If so, the scattering of the dominant "electric" particles on 
this magnetic field is an important component of transport, and this is what we are 
going to study below. 

In order to complete this brief overview, let us comment on the issue of "Hig- 
gsing" in hot gauge theories. It is well known that, at high T, one can consider 
a dimensionally reduced effective 3d theory, which is basically a Georgi-Glashow 
model with Aq being the adjoint scalar. Its lattice-derived phase diagram, as well 
as the effective lines of the high-T 4-d gauge theories derived to 1 and 2 loops can 
be found in Ref. [28j. If taken literally, the line for QCD crosses into the broken 
phase, except at extremely high T: but the authors themselves argued that this 
line should in fact remain in the symmetric phase. Thus, even the asymptotically 
high-T behavior is under debate. In the temperature range of few- Tc we know that 
the expectation value of the Polyakov line changes from zero to unit value, and 
that its effect on quarks [29] and gluons [30] is large, producing significant suppres- 
sion. This approach in particular gives a very good description of the quark number 
susceptibilities [3I|. 

One more practical aspect of the issue comes from heavy-ion phenomenology. 
Dumitru and collaborators [32J have used this form of the effective Lagrangian to 
study the real-time evolution of Aq. The main conclusion from their work is that 
{Aq) belongs to the class of so-called slow variables, and its evolution in heavy ion 
collisions has to be treated separately from the overall equilibration. They have 
numerically solved the EOM for Aq, starting from the "suddenly quenched" value 
corresponding to its vacuum form, moving toward its minimum at the deconfined 
phase at T = 2Tc. The main finding of this work is that the relaxation of this 
variable is very slow, taking approximately 40 fm/c. This time significantly exceeds 
the QGP lifetime at RHIC, which is only about 5 fm/c, which suggests that in 
real collisions we should treat Aq essentially as a random variable frozen at some 
value and color direction during hydro evolution. This means that there is a chaotic 
out-of-equilibrium Higgsing, slowly rolling down, like in cosmological inflationary 
models: thus one would like to know as much as possible about phase transitions 
and EoS for all values of the Higgs VEV. 
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2 Monopoles in gauge theories at high T 



2.1 The monopole density and properties 



In section 1.4| we have already briefly reviewed the main issues and ideas related 



with magnetic monopoles in hot QCD. Now we begin a much more quantitative 
discussion of the main parameters and lattice results involved, which will be needed 
for our discussion of transport cross sections. 

In Fig. [T] we compile the lattice-based data on the monopole density. Since 
we compare various theories, SU(2) and SU(3) pure gauge theories as well as those 
with quarks, we need to explain the units. Following the lattice tradition, physical 
units are defined by insisting that the string tension is the same in all of them, 
^/(T = 426 MeV. We will always show the temperatures in units of the corresponding 
deconfinement transition temperature Tg. 

In Fig. flla) the vacuum value of the density is taken from Bornyakov et al. [23] , 
(nm)T=o/o" = 0.5. The dots and crosses are from Liao and Shuryak [Mj: they 
correspond to "condensed" and "decondensed" monopoles, respectiveljj^ Note that 
these two add together to an approximate constant total monopole density, for all 
T <T,. 

In Fig. [l]^b) we summarize what is known about monopoles in the decon- 
fined phase. The diamonds of different colors show the direct lattice observation 
of monopole density by D'Alessandro and D'Elia pT], scaled up by the factor 2 
which accounts for the transition from SU(2) to SU(3) gauge group. The best fit to 
the data of D'Alessandro and D'Elia (not shown) is 

njT^ = A/\og{T/K^Y (13) 

with A = 0.48,0; = 1.89 and Tc/Aeg = 2.48: we discuss the expected parametric 
dependence at high T below. 

The crosses in Fig. [l](b) are from Ref. [34j, as those in Fig. [l|^a), but for tem- 
peratures larger than Tc. They have been obtained from (metastable) flux tubes 
in the plasma phase. The latter should be compared to the solid line, which rep- 
resents gluons, affected by the nonzero VEV of the Polyakov line effect^ |30] . (At 
high T the normalized gluon density becomes a constant corresponding to Stefan- 
Boltzmann ideal gas of massless SU(3) gluons.) Note that the densities of gluons 
and monopoles cross each other at T ~ 1.3Tc. 

The available lattice results on the internal monopole structure are not too de- 
tailed. Fig. |2| is taken from the talk given by Ilgenfritz at Lattice 2007 [35] . Panel 



^The origin of those points resides in the lattice results on string tensions corresponding to free 
energy and potential energy, respectively. 

^The reader is reminded that it reflects the nonzero VEV of Aq, which acts on electrically 
charged particles as an imaginary chemical potential. 
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(a) 



(b) 



Figure 2: The square of non-Abelian components of the gauge potentials (a) and 
the topological charge (b) correlated with the monopole path (taken from Ref. [53]). 



(a) and panel (b) show the structure, transversal to a monopole world line, of the 
"non-Abelianicity" (the gauge potential not in the Abelian direction) and of the 
modulus of the topological densitjj^ respectively. These two quantities are defi- 
nitely correlated with the monopole current, however they display different sizes. 
The non-Abelianicity shows a narrower peak, with a small width of only about 



1.5a ^ 0.15 fm 



(14) 



(The lattice spacing was a = 0.105 fm in this simulation). It is hardly surprising 
that we don't have any real insights about the monopole structure, as their size is 
so close to a (the UV cutoff) of the best lattice calculations. 

From Fig.(b) it is evident that the topological charge shows a larger radius: the 
interpretation is that this is actually the size of dyon^ having both E, B along the 
radius. Although we will not go into this subject, let us note that selfdual dyons, 
possessing the topological charge and zero modes, have been studied via fermionic 
methods without any gauge fixing: nice correspondence to Kraan-vanBaal solution 
[36] was found. 

The issue of monopole mass is currently under intense study. We have been 
informed by D'Elia and D'Alessandro about their derivation of the mass from 



^Provided by the unfiltered overlap-defined topological density. 

®The reader is reminded that the lattice definition of monopoles, based on the end of the Dirac 
strings, includes all magnetically charged objects. The density that we showed in Fig. [T| also 
includes all monopoles and dyons of all kinds. 
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monopole paths. We will only mention that the mass grows with T substantially, 
justifying the static monopole approximation adopted in this work. Only close to 
the monopole mass becomes lighter than the gluonic one, reaching about 300 MeV 
at Tc, as anticipated from Bose condensation arguments [5j. This region is excluded 
from discussion in the present work. 

2.2 Parametric dependence at high T 

Let us now switch to the domain of asymptotically high T and discuss what should 
happen when the electric sector gets really weakly coupled, namely when (T) /47r <^ 
1. At high T ^ Tc we are in weakly coupled QGP, with "electric" degrees of free- 
dom - quarks and gluon^ dominating. Perturbative scales include (i) the hard 
scale T; (ii) the electric scale, e.g. that of the electric mass found perturbatively [2], 
Ee ~ eT. In the magnetic sector, perturbatively there is no screening but power 
infrared divergencies of the diagrams, on which basis Polyakov [17J conjectured that 
magnetic screening will be developed nonperturbatively at the (iii) magnetic scale 

Em ~ e^T ~T/log(T). (15) 

This conjecture is supported by the lattice data at high T (e.g. by Nakamura et 
al. ^37]). To give an idea of the magnitude, let us mention that at T/T^ k, 10 the 
numerical values of the two screening masses are ]\/['^^^^^'^^'^9 jj^ ^ -^^ j^^reemng ^rp ^ 

2.3. The total contribution of the magnetic sector to thermodynamics (pressure) is 
of the order of 

Pm ~ e^N!T\ (16) 

similar to what one obtains from 4 (and more) gluonic loops. This value is much 
smaller than the total pressure p ~ N'^T^. 

Perturbative arguments for the magnetic sector suggest for its entropy, monopole 
density, magnetic screening mass and spatial string tension to scale with Polyakov's 
"magnetic scale" e^T . In fact, lattice results for the monopole density [21] that we 
discussed above can be fitted as ~ log(T)^^ for T /Tc > 4 (omitting the near-T^ 
points which are not expected to obey the asymptotic behavior). 

One general property of the magnetic scenario [2| is the expectation of a strong 
magnetic coupling at high T. More specifically, it is expected [26] to be a monopole 
plasma in liquid form, with the ratio of potential to kinetic energies (known as plasma 
parameter F) to be constant T{T) —>■ r*(oo) ^ 5. Yet, there are many more specific 

^One may ask why we call gluons "electric" excitations, as they have E and B of equal strength. 
A more accurate formulation of what we mean is: they are the excitations which can be inferred 
directly from the Lagrangian of the "electric" formulation of the gauge theory, in which the mag- 
netic objects are solitonic objects. The "magnetic" formulation would have a Lagrangian with 
monopoles as sources, containing of course the gauge fields but no quarks. 
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questions about the magnetic sector which remain open. In general, the monopole 
properties are formed under the influence of the following three important effects: 
(i) interaction with Aq, the "Higgs field"; (ii) interaction with other monopoles; (iii) 
interaction with "electric" quasiparticles, quarks and gluons. 
Let us comment on them, subsequently. 

(i) At nonzero T the role of scalar adjoint Higgs is taken by Aq, whose nonzero 
VEV breaks the color group into massive non-Abelian and massless Abelian gauge 
bosons. This phenomenon, related to the expectation of the Polyakov loop, is widely 
studied, see e.g. [3l]: it is most important at T ~ (1 — 4)Tc. The Euclidean-real 
eAo/T produces a 0(1) phase factor - also known as imaginary chemical potential 
- for electric particles, quarks and gluons. The corresponding extra i in Minkovski 
VEV does not allow to use standard 't Hooft-Polyakov monopole solution of the 
Georgi-Glashow model directly: yet one would still expect that the monopole mass 
is heavier than that for quarks/gluons by the factor l/e^(T) ~ log(T), thus with 
Mm/T growing at high T. 

(ii) Monopole-monopole interaction is Coulombic, and (as shown in [26]) the 
available lattice data on same and opposite sign monopole correlations are in good 
correspondence with MD simulations for such plasma. One important point is 
that, at high T, the stronger magnetic coupling cancels the relevant diluteness of 
monopoles, thus the corrections to their mass are finite 

AM^/T ~ /n„i/VT ~ g^e^ ~ 0(1); (17) 

similar implications hold for the magnetic screening mass 

(M™'"^)' ~ g'^njT ~ g^e^'T^ ~ [e^Tf (18) 

which is thus at the expected "magnetic scale" . Below we will consider more detailed 
models of the monopole correlations and the collective (magnetic or dual) potential 
resulting from them. 

(iii) Electric-magnetic interaction is what this paper is mostly about. Let us 
only note here that, since such interaction is based on the Lorentz force, it contains 
the coupling combination ge ~ 0(1) locked by the Dirac condition, thus here there 
is no small parameter involved. 



2.3 Parametric estimates of the scattering rates 

We define the dimensionless scattering rate of a particle of type i on all particles of 
type i as 

^ ^ ^^^^ 

where Uj is the density of "scatterers" and f jjCTjj is the relative velocity and corre- 
sponding (transport) cross section. 
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We start by counting the powers of log(T) in the well-known results for gg 
gg, ggg processes, to be compared with the monopole-related ones. 
We have 

- K^-) ^ «^ ^ f _4_y (20) 



T T ' Vog{T) 

Although the radiative process gg —>■ ggg is a higher-order process and has extra 
coupling, it also has one "soft" propagator carrying momentum at the "electric 
scale" eT, which puts it at the same level as the elastic process. 

As we argued above, the monopole density scales with the "magnetic" scale: 
nm ~ (e^T)^ and thus it is small; nevertheless, there is no couphng constant in the 
cross section. Thus 

) log(log(T)) 
T T ^ log^(T) ^ ' 

where the factor log(log(T)) comes from divergent transport cross section (to be 
derived below). Thus, the process that we consider is parametrically subleading at 
high T, but not by much. 

There are many other details to be included below - such as numerical constants 
in all densities and cross sections - but at this point we only address the dependence 
on one more general parameter, the number of colors Nc- The number of gluons 
scales as N^, while that of monopoles scales as A''^^. The perturbative gg cross section 
contains the 't Hooft coupling squared, = (e^A'c)^ with extra powers of Nc, but 
the result only makes sense as an expansion in A <^ 1, so the limit of large Nc, 
irrespective of e, cannot be taken. The non-perturbative gluon-monopole amplitude 
has partial scattering phases which are pure numbers 0(1), thus one can rather 
write that the ratio of interest is Wgm/wgg ~ 1/X^Nc, with small A but large Nc- 

2.4 The monopole mass 

In the original context of the Georgi-Glashow model, the monopole mass can be 
written as a perturbative series 



9^ 



(22) 



The first classical term contains the known functions /o(0) = l,/o(l) = 1.238, 
obtained from the numerical solution of the classical EOM. As for the second term, 
only its limit at small argument has been calculated [38] 

/i(;2) = -^log(z2) for z^Q. (23) 
/vr 

There was one lattice calculation of the monopole mass, by Rajantie [52], which 
was performed a.t z = Mh/Mw = 1, with very weak coupling = 1/5, at which the 
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deviations from the classical limit were small. From these results one can estimate 
that /i(l) ~ —(1.5 — 2). If so, extrapolating it to cig = e^/47r ~ 0.7 (which we 
have at T ~ a few times T^. in QGP), one finds that quantum corrections would 
be as large as the classical mass, cancelling it to zero. While one cannot trust this 
extrapolation quantitatively, it certainly shows the trend. 

What is the mass which is needed to reproduce the observed monopole density? 
We write the expression for the density as 

^ = Ml /■ (24) 

^"^ 27r2 J exp[(6, + V{Q))/T] - 1 

with = +p ^. We included in the above expression the mass as well as the 
correlation energy, which is due to a (dual Coulomb) potential induced by all other 
monopoles at the position of one of them. We will write it as 

v{0) ^ ( 



^ . rr. ,M (25) 

with a = (n^)^^/'^ and M the so called "Madelung" dimcnsionless constant: its 
value is the subject of the next subsection. The effective degrees of freedom are 
dof=2 for the SU(2) gauge group and 4 for the SU(3) one, which corresponds to 
(spinless) monopoles and antimonopoles. 

If one wants the monopole ensemble to mimic the scaling relations mentioned 
above, the mass (minus V{0)) should slowly grow with T 

^ 31og(log(T)) (26) 

in contrast to the perturbative masses of quarks and gluons which slowly decrease: 
rUg/T e [log(T)]~^/^. This leads to the picture that we consider here, namely 
heavy, rare monopoles embedded into a dense plasma of light gluons and quarks. 
If so, the main interaction to be included is a strong-coupling mutual magnetic 
Coulomb one. 



2.5 Monopole interactions 

Before we get into the details of the possible models describing the monopole en- 
semble, we provide the numerical values of some key variables at two temperatures 
(we choose them near the ends of the interval but, prudently, we do not take the 
endpoints) : 



T/T, 


nm(fm"'') 


aT 




^V(47raT) 


1.424 


2.834 


1.5 


0.29 


1.3 


9.865 


169.2 


2.6 


0.056 


2.4 
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Let us start by comparing the dimensionless densities n^jT'^ = 1/ ifl^^"^ iii the table, 
to that of a massless Bose gas with two degrees of freedom, n/T^ = .244 (hke that 
for the blackbody photons). While at the higher T/T^ = 9.865 the monopole density 
is several times smaller than this number (ascribed to a mass-induced suppression 
above), at the lower T/Tc = 1.424 it is somewhat above it, which would require 
"overcompensation" + ^(0) < 0. Since in this work we will not discuss the 
near-Tc region, with all its complications, we conclude for now that T/T^ = 1.424 
is a bit too low for our basic picture to apply in full. The reason why we will still 
include its discussion below is pragmatic: monopole correlations are not yet known 
at really high T. 

The value of the effective plasma parameter {g'^ /{A-KaT)) is given in the above 
table, the coupling is from correlations [21]. Note that, if the Madelung constant is 
simply=l^'', the interaction-induced factor is about exp{g'^ / (AiraT)) = 3.6, 11 at the 
two considered temperature values. This cannot be taken literally, but shows that 
the interaction effect can be quite substantial. The values of F itself grow from 1 
at T = Tc, to about 4.5 at T/Tc = 4, as calculated in ref. [26]. Its T-dependence 
shows that it approaches a constant value r*(cxo) ?a 5 for large T. 

In the expression for the density, V{0) implicitly depends on how the other 
monopoles are distributed in space, so we have to solve a manybody plasma problem. 
The distribution of those monopoles in space can be described by three models, 
correponding to (i) "weak" , (ii) "medium" and (iii) "strong" correlations. Physically, 
they correspond to a "gas", "liquid" and "molten salt" examples, respectively, all 
well known and studied for classical Coulomb plasmas. 

In the first case (i) we simply ignore the correlations (perturbative Debye theory 
includes a cloud of correlated charges with total charge -1, but with large size and 
small Debye correlation energy E = —{g"^ /{AuRd)) ~ g^T). 

The second model [IQ] (ii) is based on two assumptions/inputs: that the static 
two-body correlation between particles is given by the lattice results for monopoles 
I; and that manybody correlations are absent. If so, the correlated charge density 



is spherically symmetric and thus 

AQ(r) = {nj2) f {g++ - g+_)47Tr^dr. (27) 
Jo 

The interaction field and the potential induced by the correlated charge can be 
calculated via Gauss' law. For one typical case discussed in [40j| in detail, namely 
for T/Tc = 1-42 when am = g"^ /^t^ ~ 2 and a ~ .24 fm, the potential induced by 
the correlated charge is harmonic at small r. Its value at the origin is V^(0) ~ 35 
MeV, thus the correction induced by it is exp(— y(0)/T) ^ 1.086. Such corrections 
0(10%) are typical for this model (ii). 
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We will discuss its realistic values below. 
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The third model (iii) has two different assumptions: (a) that local multi-body 
correlations are significant; (b) that those can be described by a cubic lattice of 
alternating charges, since this is known to be the ground state for ion lattices such 
as NaCl (thus the name "molten salt"). Collective potentials in solid plasmas are 
easily calculable from obvious triple sums over the atoms. Although it is not difficult 
to calculate them numerically, as the sum over expanding cubes converges well, quite 
a lot of mathematical efforts have been made over nearly a century to get the best 
analytic representation of it. The latest one that we found is ref. [2], which tells 
that the so called Madelung constant is 

where 5* is explicitly exponentially small O(10~^°). 

Physicswise, by elevating the Coulomb energy of a particle at the exact crystal 
node to a Coulomb potential, induced by all charges but one, one can see the di- 
rections in which it decreases. This reflects the obvious fact that classically the ion 
lattice is simply unstable against opposite charges falling at each other. That does 
not happen quantum-mechanically, of course, because of the "localization potential" 
~ h'^/mr'^. Furthermore, both atoms in salts and monopoles in QGP are extended 
objects, which only interact via the Coulomb law at large distances r, while there 
exist some repulsive cores at small r. For molecular dynamics simulations [7| we 
used a potential of the form 



V{r) 



Ana 



+ 



(r/a) p{r/a)P 



(29) 



where p is some power, chosen to be 9. (The factor {1/p) in the second term is 
needed to make the force zero at distance r = a.) It leads to a nice stable in-lattice 
potential with a near-isotropic harmonic well for small deviations from the node. 
The energy per particle for this potential is 

V(0)^-1.1— (30) 

^ ^ Ana ^ ' 

So far we don't have lattice results on manyhody correlations of monopoles and thus 
we are not sure which of those is more realistic. We have presented those results 
for future comparison, and also to emphasize that the collective magnetic (dual) 
potential is large and non-negligible compared to the monopole mass. Both should 
be included in the monopole density calculations, which are thus more involved 
compared to a weakly coupled gas considered before. 
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3 Scalar fluctuations 



We set up the problem of quantum scattering on monopoles in the framework of the 
Georgi-Glashow model |12] (see Appendix A for the details on the model and the 
notations we use). The generalization to QCD will be discussed in the following. The 
Georgi-Glashow model describes the interaction between gauge fields (two massive, 
spin-1 Ws and a massless photon) and a scalar Higgs field. We start with the 
scattering of the scalar fluctuations on monopoles. We introduce a total angular 
momentum operator T, which is the sum of the orbital angular momentum L and 
the isotopic spin /: 

f=-ifxV + I (31) 

with {I°')bc = —i^abc- In terms of these operators, the wave equation for the scalar 
fluctuations can be written in the form (for a derivation of this equation see Appendix 
A): 



2d_ 

r Or 



T2_ 



f ■ / 



■X 



-dl 



A 



X + ^ —X (32) 



2^/7(0^- + 

e 



— V 



X 



0. 



The term which is proportional to induces charge-exchange reactions. 

We can define a simultaneous eigenfunction 4'u"'{'r)a of the commuting operators 



T^, T3, I and f-I (see eq. (11)). This function depends only on the angular variables 



specified by f. A solution to the equation (a specific partial wave) can be written 
as the product of the angular function 0™"(r) and a radial function S'"(r) 

Xirla = ^rir)aSr{r). (33) 



The angular function (t)u'^{r)a is peculiar because the operator I is projected along 
f. Therefore, the angular function must be rotated, from the standard cartesian 
frame, to a "radial" frame. This construction can be achieved by making use of a 
spatially dependent unitary matrix which rotates f ■ / into I-^: 



(34) 



We therefore have 



f-IU{-^,-9,^) 



(35) 
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where 

f = -f X ( i V + eALs ) + fh (36) 



and 

eA=—. (37) 

r + z ^ ^ 



Eqs. (11) are satisfied by the following function 

(l^rir)a = {U{-v,-e,^)x7)aV{f) (38) 
where {Xi)a is an eigenvector of I3 in the cartesian basis 

h{x7)a = n{x7)a (39) 

and the function V{r) obeys 

ll?(f) = ( ^(^ + 1) \v{r) (40) 



where /a in eq. (36) is now replaced by its eigenvalue, n. We have 

m = = {t,n\e-^^^^e^'^^e^^^%m). (41) 

We can write the function (f)^"'{r)a by making use of the following expansion 

X (42) 



n' 
t+i 

X ^ {i,-n,t,n\l,{])vf^,^_^,{-ip,e,ip){l,m - n'\i,-n\t,m) 
i=\t-i\ 



where 



<„.(«,A7) = V^^'"(A7)- (43) 



These functions are normalized in the following way 



Since the above function is a polynomial in x/r, y/r, z/r, the function in eq. (42) is 
analytic everywhere. This clearly shows that there are compensating singularities in 
U{—ip, —6, if) and P(f). Another way of building the spherical harmonics (42 ) would 
be to act with the vector operators r, (rV — iA), (rV + iA) on the standard angular 
functions Ytm{d, ip) (with A = r x V). The two methods are totally equivalent, but 
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the one we choose here, following ref. [TO], provides an easy way to understand the 
mechanism of rotating from the standard, cartesian frame to a "radial" one. Since 
equation (32) in its most general form admits mixing between particles of different 
charge, in order to get the equation for the radial functions S*" we have to write the 
order 1 fluctuations of the field around the classical solutions as a superposition of 
the functions describing a particle with definite charge n: 

x{r,xo)a= Yl e^""°^0r(^,¥')a (45) 

n=-l,0,l 

where the three functions 5'"(r) correspond to the three physical fluctuations with 
charge n = 0, ±1. We plug the above expansion for x(0 into Eq. (32); through 
this procedure we obtain the following system of equations for the radial functions 



tit + l) , KiO 



+ 2- 



UJ 



SKO - A 3 



1 S^iO 



+ 



(46a) 



t(t + l)-l , K{0' 



— UJ 



1 slio 



(46b) 



t(t + 1) 



e 



1 sr\o 



+ 



(46c) 



where we have introduced the dimensionless variable ^ = evr. At the end of this 
section we will discuss how to fix the scale and go to physical units. As it is evident 
from the above system, a mixing occurs in the monopole core between different 
charges: the term oc K{C,) involves a mixing between charges that differ by one 
unit. The above system of equations has been obtained in the most general case, 
for generic angular momentum t. Nevertheless, we have to keep in mind that there 
are some restrictions due to the following requirement: 

r-T = f- L + f- I = f- I = n. (47) 
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For this reason, in the case t = only the n = scalar fluctuation is allowed. The 
equation for this special case and its solution will be discussed in the following. 



For t > 0, the system of equations (46) is difficult to solve, due to the mixing 
between the different radial functions. This mixing is due to the charge-exchange 
reactions that can occur inside the monopole core. If the monopole core is small 
(we have seen in Section 2.1 that lattice-based estimates for the monopole size give 

— 0.15 fm) we can neglect the charge-exchange reactions. This corresponds to 



considering the above system of equations (46) in the limit 

K{0 ^ H{0 ^ e 

In this approximation, it reduces to 

't{t + l) 



+ 2A S°(e) = 



(4J 



(49a) 



sl"{0- 



sm = 



(49b) 



t{t 



LU 



(49c) 



From the above system it is clear that we can identify the radial functions with 
spherical Bessel functions having index t', which is the positive root of 

,2 



t'{t' + 1) = t{t+l) -n\ (50) 

Namely, in the limit of small monopole core we have S^{r) jt'{kr). In general, t' 
is not an integer number. The corresponding scattering phase will be 6t' = t'7r/2, 
independent of the energy of the incoming particle. 

For t = we have only one fluctuation allowed, namely the one having zero- 
charge: Sq{^). It obeys the following equation 



) sm - X - 1) sm = -^^^0° (0 



(51) 



In this case, there is no Coulomb potential of the form 1/^^, which is obvious since 
a charge-neutral particle does not feel the Lorentz force. The scattering in this case 
is entirely due to the monopole core. We can solve the above equation numerically, 
thus obtaining the scattering phase as a function of the energy of the incoming 
particle, by imposing the following boundary conditions 



sm 

S'o'iO 



sm 



- 2A + 5o 

^Vtu2 -2X + 60 



00 



y/uj'^ — 2X cos 



00. 



(52) 
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10 15 20 25 
1^1 

Figure 3: Scattering phase Sq as a function of \k\ = ^/u"^ — 2A. 6o is obtained by 



solving Eq. (51 ) with the boundary conditions (52 ). The continuous hne corresponds 
to the BPS hmit (A = 0), while the dashed line corresponds to A = 1. 



The scattering phase as a function of = y/u'^ — 2X is plotted in Fig. s] in the 
BPS (non-interacting) case, and also for a finite value of the coupling A in the Higgs 
potential (A = 1). Fig. |4] shows the classical solutions H(^) and K(^), both in the 
BPS limit and for A = 1. At this point we need to fix the scale in our problem, 
and to estimate the scattering length in physical units. In order to do it, we have to 
connect the physical size of the core of the monopole, (14), to the dimensionless 
units that we used in the Georgi-Glashow model. We recall that we obtained 
through the width at half height of the "non-Abelianicity" , which is defined as the 
square of the non-abelian components of the gauge potential. Therefore, it is natural 
to fix the scale by imposing that, in the Georgi-Glashow model, coincides with 
^m, the width at half height of the function K{^y. By looking at Fig. |4| we can see 
that, in the BPS limit, we have C,m = 1-49, while for A = 1 we have C,m = 0.87. Thus, 
ev = 9.94 fm~^ in the BPS limit, while ev = 5.8 fm~^ for A = 1. The scattering 
length is defined as a^z ^ —d5t'/dk. From Fig. [3| which shows the behavior of 5o as 
a function of \k\, we get agi = 1 in the BPS limit and agi = 0.6 for A = 1. Converting 
to physical units, we obtain agi ~ 0.1 fm for both values of A. 
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Figure 4: Classical solutions K{^) and H{C,) for the 't Hooft Polyakov monopole at 
A = (BPS limit) and A = 1 (taken from Ref. gS]). 

4 Vector fluctuations 

In the case of vector particles, the generalized angular momentum J is made up of 
three components: the orbital angular momentum L, the isotopic spin I and the 
spin S: 

J = L + I + S = f+S. (53) 

There are generally two different ways of composing three vectors, depending on 
the two vectors to be added first. The monopole vector spherical harmonics are 
eigenf unctions of and J3. Due to the following relation 

f ■ J = r ■ L + r ■ 1 + f ■ S = r ■ 1 + r ■ S = n + a, (54) 

the allowed values of the total angular momentum quantum number j are |^| — 
1, I n I, ... except in the case of ri = 0, where |n| — 1 is absent. There are different 
ways of building the monopole vector spherical harmonics; for example, they can be 
constructed by making use of the standard Clebsch-Gordan technique of addition of 
momenta. 

Another possibility, which we will adopt here, is to build the vector harmonics 
with j > n by applying vector operators to the scalar harmonics, as we will see in the 
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following By definition, these harmonics can be introduced as eigenfunctions 
of the operator of the radial component of the spin, S ■ r. We will show that this 
is a very useful choice for the "hedgehog" configuration we are working in: in fact, 
there is a natural separation between radial and transverse vectors, making this 
choice particularly useful for studying spherically symmetric problems. The vector 
harmonics with the minimum allowed angular momentum j = |n| — 1 cannot be 
constructed in this way and must be treated specially. As already mentioned, there 
is more than one way to obtain a given value of j, and thus several multiplets of 
harmonics with the same total angular momentum. In the following, we will classify 
the multiplets by the eigenvalue of r ■ S = a. In general, a = 0, ±1, but it is further 



restricted by the requirement (54), which implies that n + a lies in the range —j to 
j. This gives 

• for j = 0: 

— n = and cr = 

— n = 1 and cr = — 1 

— n = —1 and a = 1 

• for j = 1 all combinations are allowed, except: 

— n = —1 and cr = —1 

— n = 1 and a = 1 

We denote the vector harmonics by ^^^{9, ip)ai- They obey the following eigenvalue 
equations 



1 rj(j+i) 

^3 

(f ■ /) 



cr 



4.1 j>2 

In this case, all the possible combinations of n and a are allowed. We can there- 
fore construct a set of vector spherical harmonics ^J^^ {9,ip)ia that will obey the 
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eigenvalue equations (55) (in the following, A = r x V): 

1 



+ 



Vj(J + 1)(2j + 1) 
1 



rV + tA)faY"'{e,ip) 



rV + iA) (rV-iA) Y"'{e,ip) 



1 



vWnwTT) 
1 



rV-zA r„y™(^^,(^) 



j(j + l)v/2(2j + l) 



rV-iA) (rV + zAj Yp{e,ip) 



(2j + l)(j-l)j(j + l)(j + 2) L 
+ ^(rV + ^A) (rV + ^A) Ff(^,<^) 

^ W + iA) frV-iA) y"(^,</^) 



rV + «A f„y'"(^,(/.) 



+ 
+ 



(2j + l)(j-l)j(j + l)(j + 2) 



rV - lA) (rV - lA) Y"'{e,ifi) 



jU + i) 



rV-zA) (rV + iAj (^,<^) 
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2j + l 



v/j(j + l)(2j + l) 



Vj(j + l)(2j + l) ^ ^« ' 



(56) 



The above functions are normalized in the following way: 
Tdy. dcos(^)($^-;-(^,y.))t$;>'(^,y,) 2 



2j + l 



Ojj'Onn'O 



(57) 



We can therefore expand the vector fluctuations in terms of the above spherical 
harmonics basis: 



a{r,xo)ia 



[T.iW<^>ro'(^.V^) +^.2(r)$ri^(^,^) +T,3(r)<|.-:^,(^,^) 



+r.4(r)$7o-^(e,^) +T,5(r)$-r^(^,^) +T,6(r)$-r/(^,^) + 



m, — l/ 



; m,— 1 / 



r.7(r)$,™o°(^>¥') +^.8(r)$,7(e,^) +T,9(r)$-::^,(^,^ 



m,0, 
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The coefficients of the expansion are radial functions; by plugging eq. (58) into 



the corresponding field equation for the vector fluctuations (A-15), we obtain a 



system of differential equations for these coefficients. We find that the equations for 
Tjj{r), Tjs{r) and Tjg{r) do not contain a second derivative with respect to ^. These 
amplitudes are closely related to infinitesimal gauge transformations (see for example 
the discussion in Ref. [13]) and can in fact be removed by such transformations. 
Explicitly we have 



mo 



2 ^^riO + -f^T;,{0 



-LU 



j{j + l) + 2K{0' 



+ 



3{3 + l)K{i) 

2 e 



MO + 



2e 

j(j + i) 
2e 



V2j(j + l)^^Tj5(0 + \ T^TjriO + ^ TjsiO = 



(59a) 
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, , jU + l) -2K{if + 2H (if 

—uj -\ 



2e 



+ 



(j-l)(j + 2)ir(0 



V(j-l)j(j + l)(j + 2) 

2e 



+ 



(j-l)(j + 2) 1 

2 e 



(59b) 



2 ^ + 2 + 4X (0' + 2i/ (0'' 
—u -\ 



2e 



-K{0 



- A/2KJTI)^r,4(0 + ^^.5(0 + 



V(j-l)j(j + l)(j + 2) 



2^2 



(59c) 



j{j + l) + 2K{iy 
—UJ -\ 



2e 



^.4(0 



+ ^^T,,(o - y^7(JTI)^T,3(0 + V'^^^^-^(^) 
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r.6(0 



2 e 



(59d) 
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2e 



r,5(e) - V2j(J+i)^r,i(0 
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+ 



V(j-l)j(j + l)(j + 2) 
2^ 

2X'(0e - i^(0 



^(0 



r,-2(0 + ^f^r,.3(0 + 



3{3 + l)K{^) 



(59e) 
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(j-l)(j + 2) l^, 
2 i 



T^9 iO 



^ V(j-l)j(j + l)(j + 2) ^ 

^ ^i3l4j 



(j-l)(j + 2) X(0 

^ 



(59f) 



+ 1 
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+ 



j{j + 1) 



v/2K-(0 T,5(0 



+ 



Vj U + 1) ' V . V. . V ^ . 



(59g) 



7;'3 (0 A/2i^ (0 ^'4 (0 , /(j - 1)(2 + 3) t;, (o ^ v^ir' (o t,4(0 



(59h) 
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= (59i) 
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As we already did for the scalar channel, wc consider the above system in the 
limit of small monopole core. This corresponds to the limit 



(60) 



Besides, we gauge away the unphysical longitudinal components of our fields (Tj7(^), 
Tjsl^), TjQ{^)), after which the above system reduces to: 



2 + 



2e 



(61a) 



2e 



Tj40 + ^^^^T,^{0 = 



(61b) 



T..(0 + ^2Zl)Z|pZ+lr,,(0 = (61c) 



2 . J (j + 1) - 2 



2e 
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9 7 (?■ + 1) - 2 
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T.3(o + ^^2ZlM+M+lr^,(e).o 
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T.a(e) + ^^2^2Z|+M±lr,3(e) = (61f) 



2e 



-^'4(0 



(61g) 
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The last three equations in the above set are what remains of the equations for 
the longitudinal components of the fields, namely for Tjj{^), Tjs{^), Tjq^^), after 
we gauge these components away. They can be used as constraints for the other 
six TjjS. Indeed, after taking into account these constraints, we end up with the 
following set of (diagonal) equations for our radial functions: 

J u + 



(62a) 



(62b) 



TAO = 



(62c) 



j(j + l)-l 



^.5(0 = 



(62d) 



(62e) 



J(j + 1)-1 



(62f) 



We observe that the above equations arc those corresponding to Bessel functions 
with index / = |(— 1 + ^/(2j~^\^Yy^^^^^4n?) . In analogy with the scalar fluctuations, 
also in this case we have a constant scattering phase Sj' — /7r/2. The same set 
of equations can be obtained, as a check, if we consider the general field equations 
for the fiuct nations, in the limit of pointlike monopole core for the three fields 
Wn (n = 0,±l): 



y2 



I ■ f] +2n i S ■ f 



We remember that the operator T is given by 

T ^ L + I ^ J - S 

and therefore we have 



Wn = 0. 



(63) 



(64) 
(65) 
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In the above equation, the only nondiagonal operator is2J ■ S and we therefore have 
to understand how it acts on our angular functions. We obtain: 



2J. S<l>ll{9, ip) = v/2(j-n)(j + l + n)<l>J°(^, ^) + 2(1 + ^) 



(66) 



2J- 5$ 



if) = v/2(j + n)(j + 1 - nX'°(e, cp) + 2(1 - n)^l~\e, if). 



Since the above equations involve mixing with the longitudinal components only 



(for example ^/^ni^^ v) o^ly g^ts mixed with $j^'„(6', v^) and not with (6, (p), and 
since we have seen that the longitudinal components are spurious and can be gauged 



away, we conclude that (63 ) will produce diagonal equations for the remaining radial 



functions. In particular we get 

J2 + 5?2_2J.5 



I ■r] ±2 S -f 



j{j + 1) + s{s + 1) - 2(1 + n) -n'^ + 2n 



J^ + S^-2J -S- J-f ±2[S-r 



j{j + 1) + s{s + 1) - 2(1 -n) -n^ + 2n 



which results in the set of equations (62) for the radial functions. 
4.2 j=0 

As we have seen, for j = 0, only three combinations of n and a are allowed 

• n = and a = 

• n = 1 and o" = — 1 

• n = —1 and a = 1. 



(67) 
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rV - iA 
rV + iA 



(6^ 



We construct the corresponding set of vector spherical harmonics in the following 

way 

($[;;°(^,^)),, = V2r,rXiO,y^) 
1 

71 
1 

The above functions are normalized in the following way 

j\com{<:n{0.^))^^Z'^e,^) = 25r.^,5'"^'. (69) 

We can write the vector fluctuations for j = as a superposition of the above vector 
spherical harmonics: 

air, Xo)ia = 



(ro7(r)$°;°(^, y.) + To3(r)$°;^(^^, ^) + To,ir)^l^\9, ^))^^ . (70) 



The coefficients Toj(r) in the above expansion are radial functions which obey a 
set of coupled differential equations. The equation for the function To7(r) does not 
contain a second derivative with respect to ^. This amplitude is closely related to an 
infinitesimal gauge transformation (see for example the discussion in Ref. |15]) and 
can in fact be removed by such transformation. The remaining two radial functions 
obey the following set of equations 



To5(0 = 



(71a) 



To5(0 + 



(71b) 



(71c) 



KiO KiO 

The last equation above can be used as a constraint on To3(^) and To^{C,), in order 
to obtain two decoupled equations for these functions: 



(72a) 



To5(0 = 



(72b) 
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We can solve the above equations both in the BPS hmit and for a finite value of A 
(for example we choose A = 1 in analogy with what we did in the scalar channel), 
by imposing the following boundary conditions for ^ — > cx): 



^03(0 



sm 



1 COS 



l + 5o 



sm 



a/cj^ — 1 cos ^ Vtu^ — 1 



(73) 



We recall that the classical solutions for K{C,) and H{C,) that we plug in the equations 
are the ones that we show in Fig. |4| both in the BPS limit and for A = 1. We plot 
the potential V{^) = iH{^f + 3K{^f - l)/^^ _ ^ -^^ ^ig. |5|for these two values of A 
(we have subtracted the mass of the particle, 1 in dimensionless units); we see that 
it can in principle admit bound states. In particular, in the BPS limit V{C,) ~ 1/^ 
at large ^, namely we have a Coulomb potential which admits infinitely many bound 
states. For A = 1 the number of bound states is reduced: we still find two of them, 
as it can be seen from Fig. [7| which shows the wave function To3(^) in the case of 
negative (quantized) energy. 

By solving the above equations, we obtain the resulting scattering phase Sq as a 
function of the momentum \k\ of the incoming particle: it is shown in Fig. [ojfor the 
two values of A that we consider. As we already did for the scalar case, by calculating 
—d6o/dk at \k\ = we can obtain the value of the scattering length a^/. In the BPS 
limit, since we have infinitely many bound states, we also have an undetermined 
scattering length. For A = 1 we get d6o/dk ~ —10 in dimensionless units, which 
translated to physical units gives a scattering length asi — 1.7 fm. This rather large 
value can be explained by the presence of the two bound states: in particular, the 
second one is very weakly bound, therefore it shows a rather large effective radius. 



4.3 j = l 



In this case, the set of vector spherical harmonics coincides with (56), but two of 



them are missing, namely ^^i{6,(p) and $™l]^^(6', y?). The corresponding equations 



m, — 1 / 
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Figure 5: Behavior of V{$,), the potential in the radial equations at j — 0, in the 
BPS limit (solid line) and for A = 1 (dashed line). 



for the radial functions are in this case: 



-u^ H — ^ 



Tn(0 + ^7^13(0 



+ ^^14(0 - 2^ri5(0 = 



(75a) 



—LO -\ 



7^13(0 



(75b) 
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100 150 200 

Figure 7: Some radial wave functions in the BPS limit (left panel) and 
(right panel) corresponding to bound states in the j=0 channel. 
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1+4(0: 



^14(0 + ^^11(0 



KiOrr. , ^(0, 



-2^ri3(0 + ^ri5(0 = o 



^^5(0 



2K iCY + g (Q- 
"I 



^(0' , i^(0 



Ti5(0-2^Tn(0 



+ ^^7^13(0 + ^^14(0 = 0. 



(0 



+ (i-i^(0') 



+ i^' (0 ^ + K (0 i^' (0 ^ + ^' (0 ^ 



= 



K{0 (i-i^(0')^ + (i-i^(0') 



2^ (0 _^ ^/ ^11 (0 



+K m' iO^ + K iO' K' iO^ + K (0 K' (0 ^ 
The above system, in the hmit of large distances reduces to 



1 



Tu(0 + ^Ti4(0 = 



^I'UO-H' + i] 7^13(0 = 



ri"4(0 - 



7^14(0 + ^T^ii(0 = 
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^^5(0 - [-c^' + 1] Ti5(0 = (76d) 

TUO = (76f) 
7^^5(0 = (76g) 



Due to the constraint given by the last two equations, we find that, in the case of 
j — 1, the two radial functions 7i3(^) and 7i5(^) are actually constant. This means 
that, for this value of the angular momentum, only Tii{C,) and Ti4(^) will appear in 
the expansion of the general solution in terms of partial waves. 

5 Applications 

5.1 Density of states and thermodynamics 

Among the applications of the results obtained above, we start with thermodynamic 
quantities. Indeed, with a small admixture of monopoles in the electric plasma, one 
can use the standard virial expansion for these quantities. The Beth-Uhlenbeck- 
like formula from the theory of nonideal gases gives well known corrections to the 
partition function, induced by the scattering on monopoles 

SMm = -'^Y.(^j + l) J dk^f{k,T) (78) 
j 

where SM^ is the modification of the monopole mass, Sj{k) is the scattering phase 
and f{k, T) is the Bose (Fermi) distribution for gluons (quarks). These partial waves 
imply the decomposition of scattering amplitudes into the usual basis of angular 
polynomials 

2ikf{e, k) = ^(2j + l)(e2^^^(^) - l)P^{cose) (79) 

3 

rather than the Wigner d-functions we have to use in the current problem. One may 
think that a comphcated decomposition of the amplitude should be needed. This 
however can be avoided in general because of the following identity 
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expressing the relevant sum in terms of the amphtude itself. 

It has been shown above, that the scattering problem involves (a) many partial 
waves with j > and (b) the exceptional j = channels. In the former, the 
scattering is dominated by the Coulomb magnetic field of the monopole, while in 
the latter the Abelian magnetic field has no effect and the scattering happens on 
the monopole core. Let us discuss the j > case first, returning to the exceptional 
j = channel at the end of the section. 

Since the Coulombic magnetic field possesses no dimensional scale, one finds that 
the dimensionless combination 2ikf{6, k) has no dependence on momentum, being a 
function of the angle only. The scattering phases are just numbers: they don't have 
any logarithmic dependence on k familiar from electric Coulomb problem. Therefore, 
in this approximation the whole Beth-Uhlenbeck correction to the statistical sum 
vanishes identically! 

The origin of this somewhat unexpected result can be traced to the fact that the 
Beth-Uhlenbeck expression was derived from a semiclassical counting of the density 
of states in a large (spherical) box containing the monopole. However the semiclas- 
sical density of states, related to classical phase space, is insensitive to magnetic 
fields because the corresponding integral 

= / 95^^ ^^"^^^'""^^ ^^^^ 

for an electric particle in any magnetic field H = {p — eA[x)Y /2m does not depend 
on the field at all (in order to see that this explanation is correct, consider for 
example a scattering on a dyon, which was also solved in the same papers [TUl [12] ■ 
In this case there is a Coulomb phase 5j ~ {^iQi + ^292) ^og{k) and thus a nonzero 
change in the density of states). 

Let us now turn to the contribution of the exceptional j = channels. The 
corresponding scattering phase has a non-trivial fc-dependence, and thus a nonzero 
contribution to thermodynamics appears. 

Before we estimate the effect of this scattering on the thermal monopole mass, 
let us briefly clarify the issue of "mean gluon momentum" . When one calculates the 
typical momentum of black-body radiation, related to integrals like j d^kkf{k/T), 
one finds that those are k ~ 3T. However such estimate is not adequate for the 
integrals which appear in both the thermal mass (considered here) and the scattering 



rates (considered in the next section). Indeed, if in the expression (78) one can 
replace dSj/dk ^ asi where agi is the scattering length, assuming that kasi is a 
small parameter (to be justified in the following), then the integral is of the Id type 
j dkf{k/T) and thus the dominant k ~ IT. The same argument applies to the 
calculation of the scattering rates, since the integral is similar 

w ~ [ dnar^ [ d^kf{k/T)/k^ (82) 
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with the cross section being cr ~ (for non-exceptional channels). These dif- 

ferences between bulk and single-particle properties are basically enforced by the 
dimension. 

Let us now look at the absolute magnitude of the scattering lengths. Our results 
for A = 0, 1 for the exceptional channels show that the gluon-monopole scattering 
length is large a^^"* ~ 10a where a ~ .15 fm=.76 GeV~^ is the monopole core 
size. This happens due to shallow gluon-monopole bound states: we remind that we 
found two of them (one near-zero) for A = 1 and infinitely many in the BPS limit. 

However we don't expect to have these states in physical QCD because the 
effective Higgs selfcoupling (and thus mass) are expected to be large, as seen from 
near-Coulomb monopole correlations [21]. We already noted that currently the 
monopole structure is too closed to UV lattice cutoff to get a realistic value of 
the effective Higgs selfcoupling A. As soon as it is determined, one may solve the 
classical equations for K{S,) and H{E,) at the corresponding value of A and determine 
the scattering phase. It seems highly probable that this coupling is large enough 
to prevent formation of the bound states. One expects that, without bound states 
the scattering length is small a^™ 0.1 — 0.2 fm, analogous to that in the scalar 
channel. 

Thus indeed ka < 1 for temperatures of the order of (.19 GeV in QCD and 
.27 GeV in pure gauge SU (3) theory) and the contribution to the thermal monopole 
mass is small 



where the latter integral has only weak (logarithmic) dependence on the gluon 
mass m, provided it is small. Since this was obtained under the assumption that 
monopoles have a small size, asiT <^ 1, this contribution to the monopole mass is 
relatively small. 

At asymptotically high T this assumption gets violated and monopoles size would 
be at the magnetic scale, so that this parameter gets large ka ~ l/e^(T) ~ log(T) ^ 
1. So to say, hard thermal gluons now resolve the monopole structure. The answer 
in this regime then is the perturbative one-loop correction to the classical monopole 
mass, already mentioned in the introduction. 

5.2 Transport coefficients: the benchmarks 

In the gas approximation, all transport coefficients are inversely proportional to the 
so-called transport cross section, normally defined as 
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where 6 is the scattering angle. While the factor in brackets vanishes at small 
angles, the Rutherford singularity in the cross section, for any charged particle, 
leads to its logarithmic divergence. Since we will be comparing the gluon scattering 
on monopoles with that on gluons, let us first introduce those benchmarks, namely 
the well-known (lowest order) QCD processes, the gg and qq scatterings: 



dagg _ ^s^ + t^ + u^ _ 8 
^ ~ 367r V sH'^u^ 3tu 



dagg 9e^ (s^ + t4 + u4)(s2 + t2 + y2 



^5) 



6) 



dt 1287r sH^u'^ 

where (we remind) the electric coupling is related to as as usual: e^/47r = a^. 

While for non-identical particles the transport cross section is simply given by 
the cross section weighted by momentum transport t ~ {1 — z), for identical ones 
such as gg one needs to introduce the additional factor {l + z)/2 in order to suppress 
backward scattering as well. The integrands of the transport cross section for these 
two processes are compared in Fig. [8j The integrated transport cross sections 
themselves are given by 



= sl^ ( - 16 + 30 log ( ^ ) ) (87) 



mm 



mm 



where the smallest scattering angle can be related to the (electric) screening mass 
by 6'^j„ = 2* M]j/s. Note that the forward scattering log in the gg case has a 
coefficient which is roughly four times larger, as a consequence of the gluon color 
being roughly twice that of a fundamental quark. Note also that the gg scattering 
is significantly larger at large angles, as compared to the qq scattering. 



5.3 Transport cross section on monopoles 

We explained (already in the introduction) that the charge-monopole scattering 
is Rutherford-like at small angles: this comes from harmonics with large angular 
momenta (large impact parameters). However, in matter there is a finite density of 
monopoles, so the issue of the scattering should be reconsidered. A sketch of the 
setting, assuming strong correlation of monopoles into a crystal-like structure, is 
shown in Fig. |9} A "sphere of infiuence of one monopole" (the dotted circle) gives 
the maximal impact parameter to be used. 

As a result, the impact parameter is limited from above by some hmax-, which 
implies that only a finite number of partial waves should be included. The range of 
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Figure 8: Angular dependence of the integrand of the transport cross sections vs 
z = cos 6, where 6 is the scattering angle. The (red) solid line shows the integrand 
for the gg scattering, dcTgg/dt * s * t * {1 + z)/2/e'^ , while the (black) dashed line 
shows that for the qq scattering scaled up by a factor 4: 4 * dcrgg * dt * s * t/e^. 

partial waves to be included in the scattering amplitude can be estimated as follows 

W = (p.)n-^/V2 ~ «T ~ l/e\T) ~ log(T). (89) 

Since at asymptotically high T the monopole density rim ~ (e^T)^ is small compared 
to the density of quarks and gluons ~ T^, jmax asymptotically grows logarithmically 
with T. So, only in the academic limit T ^ oo one gets jmax oo and the usual free- 
space scattering amplitudes calculated in [TU] where all partial waves are recovered. 
However, in reality we have to recalculate the scattering, retaining only several 
lowest partial waves from the sum. As we will see, this dramatically changes the 
angular distribution, by strongly depleting scattering at small angles and enhancing 
scattering backwards. 

5.3.1 Scalar particle 

By making use of the partial-wave solutions to the wave equations that we obtained 
in Sec. |3| we can proceed to construct the scattering solutions. At large distance 
from the origin, they should behave as a plane wave plus an outgoing spherical wave. 

We consider a particle entering from z = —oo; the z component of its total angu- 
lar momentum is well defined: J3 = T3 = — (T-f) = —n. Therefore, the appropriate 
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Figure 9: A charge scattering on a 2-dimensional array of correlated monopoles 
(open points) and antimonopoles (closed points). The dotted circle indicates a 
region of impact parameters for which scattering on a single monopole is a reasonable 
approximation. 

angular functions that we must use to expand our solution are (pjl'"' {9 , if) . In the 
following, we consider only n > 0; a negative n would just produce an overall phase 
factor. We consider the following partial-wave sum 

Jmax 

j=\n\ 

= [u(-^,e,^)x7]^^^\f1- (90) 

We recall that the index j' in the above formula is the positive root of 

j'{f + l)^j{j + l)-n'. (91) 

The square brackets in the above formula contain a vector which is the charge 
eigenstate. Since in our approximation we don't have any charge-flip reactions (they 
would be possible inside the monopole core, which we consider pointlike), this factor 
will remain the same for the incoming and outgoing waves. The physically significant 
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wave function is therefore given by 



j=\n\ 

Jmax 



= e-^™ Yl + l)e'^'e-'^^'/%ikr)e-''-^dli;U9). (92) 
i=l"l 

The above expression contains the functions which can be written via the 

generahzed Rodriguez formula as {z = cos 9) 

So, the first term with j = n is simply ~ (1 — cos^^)", a function strongly peaked 
backwards at large n. 

Obviously, a tendency to scatter backwards is very important for transport prop- 
erties and equilibration of electric plasma. 

The asymptotic value of the wave function il)'^^\r) is then 



with the scattering amplitude f{0) given by 



pikr 

e'^' + f{e) 

r 



(94) 



Jmax 



2ikf{e) = J2 (2J + l)e^-(^'-^)rf^;L(^). (95) 
j=H 

The above sum would be badly divergent if we would sum over j up to cxo. It is 
possible to improve its convergence by separating the most singular behavior (see for 
example the discussion in [ID])- We do not have such convergence problems, since 
we have an ultraviolet cutoff on j due to the finite density of monopoles. Therefore 
we calculate the scattering amplitude by directly using Eq. ( [95| ). 

The integrands of the transport cross section (1 — cos ^^) 1/(6*) p are shown in Fig. 



To] for n = 0, j.^ax = 2,4,6 (left panel), n = ±1, j^nax = 2,4,6 (right panel). 
One can see how much their angular distribution is distorted. Strong oscillations of 
this function occur because we use a sharp cutoff for the higher harmonics, which 
represents diffraction of a sharp edge. This edge in reality does not exist and can be 
removed by any smooth edge prescription (for example a gaussian weight). However, 
we further found that the transport cross section itself is rather insensitive to these 
oscillations, and thus there is no need in smoothening the scattering amplitude. 



The transport cross section as a function of jmax is shown in Fig. 11 it is large and 
smoothly rising with the cutoff. 
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Figure 10: Scaled Integrand of the trasport cross section g(ff) = k'^{l — cos6')|/(^^)p 
with only 2, 4 and 6 lowest partial waves included, for a scalar particle with n = 
(left) and n = ±1 (right). The curves can be easily recognized by higher j^ax having 
more oscillations. 
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Figure 11: Normalized transport cross section as a function of cutoff in maximal 
harmonics retained, for n = 0, 1. 
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5.3.2 Vector particle 

We consider a vector particle entering from z — —oo. As we already saw in the 
case of the scalar particle, the z component of its total angular momentum is well- 
defined: we have seen that, if we consider angular functions that are simultaneous 
eigenfunctions of (T • f) and {S ■ f), the equations for the corresponding radial 
functions are diagonal. This means that, once we neglect the monopole core, we 
can consider an incoming particle with definite charge and radial polarization, and 
these quantities will be conserved in the final state, after the scattering. The third 
component of the angular momentum will then be 



J8 = - 



(f-f^ + (^S-f^ = - [n + (t] = -I/. (96) 

We have seen that, after gauge fixing, for each particle of given charge we have two 
possible spin polarizations a = ±1. This means that, for a vector particle with 
charge n — 1(— 1), we can have u — 0,2(0,-2) respectively, while for the vector 
particle with n — (the photon) we have u — ±1. We consider the following 
partial- wave sum 



J' 



max 



= [s(^,(^)::]*w(f) (97) 

where the factor [S (^, c/?)^] is a two-index tensor representing the charge and {S ■ r) 
eigenstate, which is not altered in the peripheral collision that we are considering. 
As we already did for the scalar channel, we consider the physically significant wave 
function 

Jmax 

*(+)(f) = e--^5](2j + l)e-^e--^'/2^,v(fcr)PgL(-^,e,(^) 

Jmax 



for which we can write the asymptotic behavior 



^(+)(^ ~ e 



(99) 



with the scattering amplitude /(^) given by 



2z^/(^)„,, = {2j + l)e'''^^'-^Ui%{e). (100) 
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Actually, there are some exceptions to the above expression for f{9): in Sec. |4]we 
have seen that the cases j = and j = 1 are somehow exceptional. For j = we 
have only two modes, with charge ±1. The corresponding radial function obeys a 
differential equation which is not Bessel-like.Yet, this solution still has a sinusoidal 
behavior at large distances, with a momentum-dependent scattering phase. We 
have also seen that, for j = 1, the only surviving modes are the two transverse 
polarizations of the particle with charge n = 0. Therefore, the above general formula 



for the scattering amplitude (100) must be separately written, for particles with 



different charge. For the neutral particle (with z/ = ±1) we have 



Jmax 



2^A;/(^^)o,±i = 5^(2j + l)4^Ii(^), (101) 
i=i 

where the factor e^'^'^^'~^^ is absent in this case, since for a neutral particle we have 
i' = j. We have seen that, for the positive (negative) charge particle, we have two 
possible values of z/ = 0, 2 (—2). Therefore we have: 

Jmax 

2ifc/(^)±i,±2 = 5^(2j + l)e^-(^'-^)cigl2(^), (102) 

3=2 



where 5o{k) is the scattering phase plotted in Fig. |6j The integrands of the transport 
cross section (1 — cos{9))\f{9)n,u\'^ are shown in Fig. 12 for jmax = 6, n = 0, ±1 and 



the corresponding values of u. One can see how much their angular distribution is 
distorted. The resulting transport cross section, 

{at)n,u = j dcose{l-cose)\f{9)n,u\^ (103) 



is shown in Fig. 13 as a function of jmax for the different possible combinations of 
charge and polarization. As we can notice, it is large and smoothly rising with the 
cutoff, in spite of the large oscillations in the integrand, due to our choice of a sharp 
cutoff. In the following we will use jmax = 6. Now we proceed to evaluate the 
scattering rate of gluons on monopoles. We recall that it is defined as 

_ {^mj'^t) gm) (104) 

where the (...) indicates an average over the incoming gluon. Following Ref. [30], we 
take into account the interaction of transverse gluons with a temporal background 
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Figure 12: Integrand of the transport cross section g(9) = (1 — cos(9))\f{9)\^ with 
only 6 lowest partial waves included, for a vector particle with n = 0, v = ±1, 
n = ±1, z/ = and n = ±1, u = ±2. 
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Figure 13: Normalized transport cross section as a function of cutoff in maximal 
harmonics retained, for n = 0, ±1 and the corresponding polarizations. 
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gauge field which is related to the Polyakov loop through 



-tr 



Pexp ( i / Aodr 



-tr exp 



iA^^Xa 

T 



(105) 



(the above definition of the gauge field actually includes the electric coupling e). If 
we select the abelian gauge for Aq, for example along A3, we get: 



A^ 

2 cos 9:^ + 1 
T 



(106) 



The presence of this background gauge field breaks color symmetry. The eight gluons 
interact with in different ways, depending on their color charge. For example, 
since we chose ^0 along the abelian direction, the two abelian gluons don't interact 
with ^0 at all. In particular, the gluon density has the following form 



8tt 
(M 



k'^dk 



+ 



+ 



exp (Pek) - 1 exp (/Je^) exp {iPA^) - 1 

1 

+ 



exp {(3ek) exp {—i/3AQ) — 1 exp {/3ek) exp (2z/3^[]) — 1 

k'^dkpg{k, T) 



exp {(3ek) exp {-2i(3Al) - 1 J i'^'^ f 



4:71 



(107) 



where we use for the expectation value which fits the available lattice results for 



the Polyakov loop in pure-gauge SU{3) gauge theory (see Fig. 15) and we define 
efc = y^/c^ + mg{TY, with an effective gluon mass rng{T) = e{T)T / \/2. For as{T) 
we use the estimate given in Ref. [26] . Notice that pg contains a factor 2 to account 
for gluon spin polarization. The coupling of gluons to the Polyakov loop gives a 
suppression of these degrees of freedom approaching from above. This is evident 
from Fig. [l| where we plot ng{T)/T^ as defined in eq. (107) as a function of the 
temperature (continuous line). 

In order to obtain the formula for the (scaled) scattering rate, we need to average 
(Tt{k) over the momentum of the incoming gluon: we use the distribution taken from 
eq. (107); in taking the integral over momenta, we have to recall that the transport 
cross section, and the gluon screening due to ^0, both depend on the charge of 
gluons: for this reason, we have to split Pg{k,T) into a "charge-neutral" part, which 
will be used as a weight for {crt{k))o,i, and a "charged part", which will be used to 
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weight {(Tt{k))ifl and (crt(fc))i,2: 



gm 



T 



+ 



+ 



+ 



n,{T)T (27r)3 J 
2 



^{^t{k)\ 
exp {(3ek) - 



exp (/5efc) exp {-il3Al) ' 
1 

exp (/3efc) exp {-2ij3Al) 



exp (/Se/c) exp {ijSA^ 
1 

- + 



exp (/?efc) exp (2i/?^|]) — 1 



(108) 



We recall that the momentum dependence of the transport cross section is trivially 
~ l//c^, except for the case (cr*(fc))i,o for which we have the exceptional case j = 0, 
with a momentum-dependent scattering phase 5o(fc). In the above equation, rimiT) 
is the monopole density as a function of the temperature. We take this information 
from the available lattice results for this quantity, plotted in Fig. [1} We show 
/T in Fig. 14 (the red, continuous line in the left panel). Also shown is the 



gm 



same quantity for the gg scattering process (black, dotted line), obtained through 
the following equation 



99 

T 



1 
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Ankldki f 2Tikldk2 



f27r)3 



-1 



d cos 0crlg{ki, ^2, cos 9)pg{ki, T)pg{k2, T) 



where agg{ki, k^, cos i 



was defined in eq. (87), with 
s = 2^1^2(1 - cos^) + 2mg{Tf 



(109) 
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6 Conclusions and discussion 

In this paper we studied the role of magnetic monopoles in the "electric" QGP 
phase. In distinction with papers by Liao and Shuryak, we have not focused on the 
near-Tc region, in which monopoles seem to be dominant over gluons in number, 
and may even expel electric fields into flux tubes as they do in the confined phase. 
Instead, we considered the more traditional QGP away from T^, dominated by the 
usual "electric" constituents - quarks and gluons, the magnetic monopoles being 
subleading in number. 

More specifically, we focused on scalar-monopole and gluon-monopole scattering. 
Classical molecular dynamics, studied in refs. [7] , suggested a very interesting mech- 
anism of "mutual trapping" by electric/magnetic quasiparticles. Based on results 
of quantum scattering on monopoles (TUl [121 US] for scalar-monopole and spinor- 
monopole scatterings, in this paper we studied this effect further. Indeed, we found 
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significant large angle and even backward scattering for small impact parameters, 
complemented with a Rutherford-like forward scattering peak. 

Our main result is a detailed derivation of the gluon-monopole scattering am- 
plitude, especially in the lowest partial waves. For gauge particles, the existence of 
spin and isospin leads to complications which so far precluded this problem from 
being solved. 

Our first important finding is that the gluon-monopole scattering has little effect 
on thermodynamics, because scattering phases - albeit large - are mostly energy 
independent, except the "exceptional channel" with j = 0. This explains why lattice 
studies which focused on thermodynamic observables were unable to see such effects. 

Our second (and the main) finding is that the contribution of gluon-monopole 
scattering is very important for transport properties. While the monopole density 
may be small, the gm scattering amplitudes have e'^g^ ~ 0{1) coupling instead of 
small -C 1. Furthermore, in our setting (with a limited number of partial waves 
j < jmax included) there is an additional enhancement for large angle (or even 
backward) scattering. 

Convoluting the cross sections found with the monopole density and gluon mo- 



mentum distribution, we plot the scattering rates nat per gluon vs T in Fig, 14 

It follows from this comparison of the gluon-monopole curve with the gluon-gluon 
one that the former remains the leading effect till very high T, although asymptot- 
ically it is expected to get subleading. This maximal T expected at LHC does not 
exceed 4Tc, where rj/ s ~ .2. This value is well in the region which would ensure hy- 
drodynamical radial and elliptic flows, although deviations from ideal hydro would 
be larger than at RHIC (and measurable!). 

The approximate relation of these rates to viscosity /entropy ratio is 

(111) 

s 5w 



We plot ?7/s in the right panel of Fig. 14 



Let us now comment on some details of the calculations, which were not yet taken 
into account and should be included later. One important fact is that the lattice data 
for the monopole density we used actually include all magnetically charged particles, 
including dyons . Their electric charge would lead to additional Rutherford-like 
scattering which we have ignored. 

In the calculation above we have not included Bose enhancement for scattered 
gluons. This effect adds the factor {l + f{p')) for gm and {1 + f{p[)){l + f{p2)) for gg, 
with the prime marking the secondary gluons. If the gluon mass is small (high T) 
those corrections are small: their magnitude is (/) ~ (T/m) ~ l/e(T) 3> 1. In the 
experimentally relevant region, when m/T = 0(1), the effect is not enhanced and 
is additionally suppressed by the expectation values of the Polyakov lineg^~ {L)'^- 
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We do not agree with Hidaka and Pisarski |47j in their conclusion that (L) < 1 makes viscosity 
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As we have already mentioned in the introduction, Xu, C. Greiner and Stocker 
[8] have suggested an alternative explanation for small QGP viscosity, namely the 
next-order radiative processes, gg ggg. Using perturbative matrix elements and 
Us = 0.3. .0.6, they found 7]/s several times smaller than for the gg ^ gg process, 
close to what we get from the gm scattering. Obviously, both mechanisms, albeit 
having such different origin, would thus be sufficient to explain the well-known 
hydrodynamic results for radial and elliptic flow at RHIC. 

It will require much more work to see how both results will change, when further 
refinements are performed. We have discussed those for monopoles above: let us 
now mention a few questions for gg ^ ggg : 

(i) Xu et al used near-massless perturbative gluons: while in RHIC-LHC range the 
lattice quasiparticle masses are instead much larger than T, about 3T or so. This 
would suppress emission of extra gluons. 

(ii) in RHIC-LHC range one should include the suppression by the Polyakov VEV 
(L) for any gluon effects (see Fig. [l] 

(iii) Inclusion of higher order corrections in badly divergent perturbative series needs 
further studies. As shown years ago in [IH], similarly treated processes gg ng with 
larger n = 4, ... lead to even larger rates! The development of convergent series for 
T]/s itself still remains to be an open challanging problem. 
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Appendix 

A The Georgi-Glashow model 

The Lagrangian of the Georgi-Glashow model describes coupled gauge and Higgs 
fields and reads 

d = -^TrF^.F^^ + Tr(Z}»(D^0)- 

= -\f;,F'^^'^ + \ (D^r) {D,r) -V{<P), (A-l) 

where F^^ = F^J"", = and Tr(J°J^) = \5ab, a,b = 1, 2, 3. For the gauge 
group SU{2), the matrices la satisfy the algebra 

[r,/"] =teatcF. (A-2) 

The covariant derivative is defined as 

D^ = d^ + leA^ (A-3) 

which gives, for the Higgs field, 

D^cp = d^<P + ie[A^,<i>] or D^cf)'' = d^^ - eea^cAlcf)' (A-4) 
and the potential for the scalar field reads 

V{ct>) = -^{ct>^r-v'Y. (A-5) 

In the above formulas, e and A are the gauge and scalar coupling constants, respec- 
tively. The field strength tensor is 

F;, = d,Al - d^Al - eeatcA^Al (A-6) 

or, in matrix form, 

Ff,, = d^A, - d,A^ + le [A^, A,] = - [D^, D,] . (A-7) 
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The field equations corresponding to the Lagrangian (A-l) can be written as 



j^^-pa^. = -ee,,e0'/^^0^ D.D'^r = -H" (^V - v^) . (A-8) 

One looks for a classical solution that is invariant under a combined space and 
isospin rotation generated by the J spin operator J = T + S, with T = L + I, and 
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S is the spin of the particle. An ansatz with this generahzed "spherical symmetry" 

is Oils] 



~amn r. 



[1-K (0] , = 0, (A-9) 

fer. These 



were H{^) and -^'(0 c'-^'g functions of the dimensionless variable ^ 
functions are solutions of the following system of coupled equations 



,(fK 



de 



KH^ + K{K^ - 1) 



,cPH 



X 



2K^H + -^H [H^ - (A-10) 



with boundary conditions 



1, 

0, 



as 
as 



^ — > oo. 



(A-11) 

In the Bogomol'nyi-Prasad-Sommerfield (BPS) limit, namely for A = 0, the equa- 



tions (A-10) have an analytic solution in terms of elementary functions 



^(0 



sinh^ 



(A-12) 



We are interested in fluctuations around the classical solutions (A-9) in the general 



case in which A 7^ 0. We therefore write the fields in the Lagrangian (A-1) as the 



sum of the classical solution, considered to be of order 0, and a fluctuation field 
which is considered to be of order 1: 



(A-13) 



The Lagrangian is then expanded up to order 2 and is therefore quadratic in the 
fluctuations and x": 

C^^'^ = -\{d^ald^al-d^ald,al-d,ald^al + d,ald,al 

-\- e \AA^Aj^(i^(iy — lA^AyCL^CLi, + 2A^A^(Xj^(iy — 2A^Aj^cb^aj^ — 2A^A^a^^aJ^ | 



- \ [2x"x" [^^^^ - t;^) + 4 ($ V)^] (A-14) 
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so that their equations of motion are hnear. 

If we only allow fluctuations in the scalar and vector channel separately, we get 
the following equations 



d^d^al - d^dX^, + ee,,, [2 (S^a^ + {d,al) A'^ + 1a\d^ A'^ - A^ [d^al) - a^d^A'J^ 
[2AlAlal - A^A^al + A^Ala^ - A^Alal - AlA^a^ - al^'^' + aj$^$'^] = 0. 

(A-15) 

The equation for the scalar particles can be rewritten in polar coordinates by making 
use of the following definition 



-ir X V; 



and by remembering that 



_2 d _L' 



dl 



(A-16) 



(A-17) 



Therefore, we can write 
92 2d 
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- A [2$''$V + x'' =0. (A-18) 



Now we remember the classical solution for the monopole ( A-9 ) so that we can write 

. ^ackLk 



-2eta,cA'^{d,xl 
[AlA^x' - A\A\x^\ 



-2i- 



Sab . r'^r^ 



and by recalling that the isospin operator is {I""] 



+ 



be 



[1-^(0]' 



(A-19) 



-ieabc we get 



_2d _ [L + I[l-K{d] 



x+ 5 —X 



-A 



g2^4 



g2 ry'2i 



(A-20) 
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The above equation can be rewritten as 

Qj,2 J. Qj. J.2 



x + 



2K{i) (l-f-{f-l) 



my 



-A 



— V 



0. 



(A-21) 



The term oc K{^) in the above equation induces charge-exchange reactions. The 
term proportional to K[^Y, modifies the deep-scattering amphtude for the charge- 
preserving processes. 



B SU(3) 't Hooft-Polyakov monopoles 

Although in the text above we deal with monopoles in a simplified way, we briefiy 
review the elements of the SU{?>) Lie algebra and monopoles. It is given by a set of 
traceless Hermitian 3x3 matrices 

T" = A72 a = 1,2.. .8 (B-1) 

where A° are the standard Gell-Mann matrices. They are normalized as TrT^-T^ = 
6"'^/2. The structure constants of the Lie algebra are f"''"^ = |Tr [A", A^JA'^ and in 
the adjoint representation we have {T"-)bc = f^^^. We are interested in the diagonal, 
or Cartan, subalgebra of SU{?>). It is given by two generators: 

/100\ /100\ 
i/i = T3 = - -1 h, = T^^—\qIq\ (B-2) 

^\000/ ^^3^00 -2/ 

— * 

which are composed into the vector H = {Hi, H2). The number of positive roots is 
3, since the dimension of the group is 8. We can take the basis of simple roots as 

/3i = (l,0), /32 = (1/2,-^/2). (B-3) 

The third positive root is given by the composition of the first two roots /^a = 
Pi — P2 = (1/2; \/3/2). Since all these roots have unit length, our choice corresponds 
to the self-dual basis P* = Pi. 



55 



Let us now remind the usual construction of a 't Hooft-Polyakov monopole in 
QCD, once we assume that the role of the Higgs is played by the expectation value 
of the field Aq, here considered to be real in Minkowski time. The QCD Lagrangian 

^= -\ [A] + d;Vuc - Di'Q,, - gr'v'.Qt] 

X [Fr [A] + Dii^Q''^ - D:,Q^'' - gfadeQK] - V[A^l (B-4) 

where we have written the gauge field as the sum of a background field and fluctu- 
ations 

K = A^ + Ql 

A^j^ background field 

fluctuations (B-5) 

and the potential V^[^o] will be defined later. In our case, the background field will 
have the following form: 

= 6,oS^'Al. (B-6) 

The field Aq will play the role of a Higgs field. Since we suppose that our background 
field is diagonal in color space, we can rewrite it as follows: 

Aq = Alh-H (B-7) 



where H was defined after Eq. (B-2) and clearly h = (1,0). If the monopole so- 
lution obeys the Bogomol'nyi equations, in the direction chosen to define the 
asymptotic magnetic field of a BPS monopole is of the form 

5„ = ^-#^. (B-8) 

Here, the magnetic charge g = g ■ H is defined as a vector in the root space. The 
charge quantization condition may be obtained from the requirement of topological 
stability, namely the phase factor must be restricted by 



exp ^ieg ■ H j = 1. (B-9) 
A general solution of this equation is given by 

9 = ^Y.''^^* = Y i""'^' + ""'^0 = + ^'^2 (B-10) 

i=l 

where rii and n2 are non- negative integers and gi, g2 are the magnetic charges 
associated with the corresponding simple roots. Our basis of simple roots is self- 
dual: PI = Pi] Pl= P2. Thus we get: 



g = g. H = — 

e 



{ni,^)Hi-^n2H2 



(B-11) 
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The magnetic charge is not so trivially quantized as in the SU{2) model. 

Since our Higgs vector h is not orthogonal to any of the simple roots jSi (B- 
[s), there is a unique set of simple roots with positive inner product with h, Thus, 
the symmetry is maximally broken on the maximal abelian torus U{1) x f/(l). In 
this case, both the numbers rii and n2 have the meaning of topological charges. 
The topological charge of a non- Abelian monopole is given by the integral over the 
surface of sphere 5*^ 

G = ^J dSnTr {BnAo) =g-h = gi + ^g^. (B-12) 

Therefore, if h is orthogonal to one of the roots f3i, only one component of g may 
be associated with the topological charge. In our case, there are two topological 
integers that are associated with a monopole. 

The above definition of topological charge can be used in order to generalize the 
Bogomol'nyi bound for the SU{N) monopoles. If we do not consider degrees of 
freedom that are related with electric charges of the configuration, it becomes 



M = Al\G\ = ^J2''^{h- a) = J2n^M, = ^ + y) 

2=1 1=1 



where = and we suppose that the orientation of the Higgs field uniquely 

determines a set of simple roots that satisfies the condition h- f]i > for all i. Thus, 



it looks like there are r individual monopoles of masses Mj. In Figure 15 we show 
the ^0 that we obtain through a fit of the lattice results for the Polyakov loop, as 
a function of the temperature. 
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